Quantum critical properties of the Bose-Fermi Kondo model in a large- N limit 
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Studies of non-Fermi liquid properties in heavy fermions have led to the current interest in the 
Bose-Fermi Kondo model. Here we use a dynamical large- N approach to analyze an SU(JV) x SU(kJV) 
generalization of the model. We establish the existence in this limit of an unstable fixed point when 
the bosonic bath has a sub-ohmic spectrum (|w| 1_e sgncj, with < e < 1). At the quantum critical 
point, the Kondo scale vanishes and the local spin susceptibility (which is finite on the Kondo side 
for k < 1) diverges. We also find an u/T scaling for an extended range (15 decades) of u)/T. This 
scaling violates (for e > 1/2) the expectation of a naive mapping to certain classical models in an 
extra dimension; it reflects the inherent quantum nature of the critical point. 

PACS numbers: 75.20.Hr, 71.10.Hf, 71.27. +a, 71.28.+d 



Heavy fermion systems show unusual quantum criti- 
cal properties, which have motivated the interest in how 
the Kondo effect is destroyed near a magnetic quantum 
phase transition |lj, LJ Li LI One theoretical approach 
studies Kondo lattice systems through a self-consistent 
Bose-Fermi Kondo model[l|. Here, a local magnetic mo- 
ment not only interacts antifcrromagnctically with the 
spins of a conduction electron bath, it also is coupled to 
a dissipative bosonic bath. (For work in a related context, 
see Ref.Li) The conduction electrons alone would lead to 
a Kondo singlet for the ground state and spin-^ chargc- 
e Kondo resonances in the excitation spectrum 6] . The 
bosonic bath - characterizing the fluctuating magnetic 
field provided by spin fluctuations - competes with the 
Kondo effect; a sufficiently strong bosonic coupling de- 
stroys the Kondo effect. A quantum critical point (QCP) 
arises, where the electronic excitations are of non- Fermi 
liquid nature. When such a Kondo-dcstroying criticality 
is embedded into the criticality of a magnetic ordering, 
the critical theory becomes distinct from the standard 
paramagnon theory |l|. This picture has been called lo- 
cally quantum critical. Direct evidence for the destruc- 
tion of the Kondo effect has recently emerged in the Hall 
measurements of a heavy fermion metal |jj. 

Previous studies of the Bose-Fermi Kondo model are 
based on an e -exp ansion renormalization group (RG) 
analvsis[lL Isl Ifj. lift 111! Il2| [the definition of e appears in 
Eq. J2J) below; e > corresponds to sub-ohmic bosonic 
bath]. This approach has been successful, partly due to 
the fact that the linear in e contribution to the anomalous 
dimension of the local s pin operator turns out to be exact 
to all orders in e |lClllllll2| . In spite of this, the approach 
is after all perturbative in e; moreover, it is capable of 
calculating only a limited number of physical quantities. 
(The model lacks conformal invariance, making boundary 
conformal field theory inapplicable.) For the case with 
an Ising spin-anisotropy, a self-consistent version of the 
model has been studied[Kj| by a Quantum Monte-Carlo 
method|14j. With spin rotational invariance, the only 
approach that has been used, beyond the e-expansion, 



is the condensed slave-boson mean field theory |l5|: here, 
the QCP is trivial: the effect of the Kondo interaction 
disappears completely as soon as the static slave-boson 
amplitude vanishes. A dynamical method is necessary to 
capture the quantum critical properties. 

In this paper, we study the model using a dynami- 
cal large- N mcthod|16l Il7l ll8J . The spin symmetry is 
taken to be SU(iV) and the number of conduction elec- 
tron channels M = kN. The large-iV limit is expected 
to capture the quantum-critical behavior of the physical 
cases (see below). This limit, with purely Kondo coup ling: 
to fermions was extensively studied in Refs. |l6lll 7l[ an d 
with purely bosonic coupling in Ref. [ig| (see also |fj,|lfj|). 

Our primary results are three-fold. First, we establish 
the existence of a non-trivial QCP even in the large-iV 
limit. Second, we show that the QCP is non-Gaussian 
not only for small e but also for e > 1/2. This demon- 
strates the inherently quantum nature of the QCP: stan- 
dard mapping of a quantum critical point to a classical 
critical point at extra dimensions |23j would imply that 
e = 1/2 is the upper critical dimension ( c-f- the classical 
spin chains with l/r 2 ~ e interactions[19|, |20|); the differ- 
ence arises from the Berry phase of the quantum spin. 
Third, we determine the critical exponents and ampli- 
tudes and the universal scaling functions near the QCP, 
as well as the properties on the approach towards the 
bosonic fixed point. We will restrict to n = M/N < 1 for 
which the spin susceptibility of the multi-channel Kondo 
fixed point is finite [m HI l2ll]. 

The Hamiltonian for the model is 

Wmbfk = (Jk/N) 2_^ S • s a + 2_^ E p c paa c paa 

+ (. 9 /Viv)s-# + ^ Wp #;-* p . (i) 



Here, a local moment 
and bosons & 

= 1 , . . . , N and a 



S interacts with fermions c paa 
p . The spin and channel indices are 



a 



= 1, 



M, respectively, and 



* — E P (*p + *- P ) contains N 2 — 1 components. We 
will first consider a flat conduction electron density of 



states and a subohmic bosonic spectrum. Denoting Qq = 
-(T r c CTQ (r)4 Q (0)) and £$ = (T T $(r)$t( )) for each 
component, we have N (u) = — jlTnQ (uj + i0 + ) = N a 
(\lo\ < D = 1/2N ) and A^(cu) = ±TmGg,(u + i0+) = 
J]p[<5(w - w p ) - d(u> + w p )} being 



h = 



A»(w) = [Al/r(2 - eJlM^sgnM, 



(2) 



for |w| < A = 1/to (r is the gamma function). Wc con- 
sider the conduction electrons in the fundamental rep- 
resentation of the SU(7V)xSU(M) group, and the local 
moment in an antisymmetric representation whose Young 
tableaux is a single column of N/2 boxes. We can then 
write S in terms of pseudo-fermions, S aa i = fj.f a ' — 
6 a ,a'/2, with the constraint J2a=i fife = N/2, which 
is enforced by a Lagrange multiplier i/z. The Kondo cou- 
pling, (Jif /iV) E CTCT ' {flier' -<W/ 2 ) C L' C ^> for each 
channel a, will be decoupled into a B^ J^ a C L/^/ 



TV in- 
teraction using a dynamical Hubbard-Stratonovich field 
B a (r) . It is convenient to consider g and Jk as indepen- 
dently varying, since the universal properties only depend 
on the ratio g/T K . 

Renormalization group analysis: Wc will first establish 
that a non-trivial QCP survives the large- N limit. Wc 
have generalized the RG equations of Ref . [ll] to arbitrary 
N and M. In the large-iV limit, the RG beta functions, 
perturbative only in e, become 



P(g) = -2g (e/2 -g + g 2 - K f 



(3) 



where, g = (Kog) 2 , and j = NqJk- There is an unstable 
fixed point at {$*,]*) = (e/2 + (1 - K)e 2 /4,e/2). Here, 
x(r) ~ 1/|t|'' and the anomalous dimension r/ = e. 

Saddle-point equations: Our primary objective is to 
study the saddle-point equations of the large- N limit, 
which have the following form (c/. Fig.QJ: 



3 W 



G 



G f (iuj n ) 



2/(r) 



1/Jk - Ss(«o;„); Sb(t) = —Ga{T)Gf{— t) 

icu n - A - T,f(ioj n ); 

k8o(t)G b (t) + 5 2 G/(t)^(t), (4) 



together with a constraint equation, 

G f (r = 0-) = (1//3) £ G / (iw„)e <w - 0+ = 1/2. (5) 



The dynamical spin susceptibility is 

x ( T )=~G f (r)G f (-r). 



(6) 



Here, A is the saddle point value of i/i, Gb{t) = 
(T T B a (r)Bt( )) and G,(r) = -<T T / [T (r)/t(o)). 

The first and second terms on the RHS of the last 
line of Eq. J3J reflect the Kondo and bosonic couplings, 
respectively. They capture the competition between the 



FIG. 1: The large- iV Feynman diagrams for the self-energies. 



two types of interactions. To proceed, we seek for scale- 
invariant solutions of the following form at T — 0, 

G f (r) = ^/Irl^sgnr, G B (r) = B 1 /\rf\ (7) 



for the asymptotic regime r 3> 
saddle-point equations become 



tq. Using these, the 



Gf(u>) 



Gg\u>) 



w + kNoB^^uj 01 



AiC ai -i 
-(K g) 2 A 1 C 1 -e +ai e> 1 - e+ai 



(8) 



,l-/3i 



Bi/9iC ft 



= ^--S B (0) + iVoAi(ai + l)C ai+ iw ai (9) 

•JK 

where C , = _»_£2M2z£M w p k avp nspf i \ . 
wneie ^-i — r(:r) sin [ 7r ( 2 _ x )/2] ■ vvc ndvt ubca A + 

E/(0) = 0, which follows from Eq. (j^J. We have also 

assumed < ai,/3i < 1, which turns out to be valid in 

most cases; exceptions will be specified below. 

Three solutions arise depending on the competition be- 
tween the Kondo and bosonic coupling terms. A dom- 
inating Kondo or bosonic term leads to the Kondo or 
bosonic phase, respectively. When the two terms are of 
the same order, the critical fixed point emerges. 

Multichannel Kondo and bosonic fixed points: When 
[3\ < 1 — e + Q'i, the Kondo coupling dominates over the 
bosonic coupling on the RHS of Eq. JHJ. The leading 
order solution is then identical to that of the model with 
only a Kondo coupling[16|. The leading exponents are 
a\ = 1/(1 + k) and /?i = k/(1 + k). 

In the opposite case, with fti > 1 — e + cti, the bosonic 
term dominates over the Kondo term. Matching the dy- 
namical parts of Eq. JHJ leads to 



ai = e/2; A 2 



2 — e ire , 

— ; -rrtan— . (10) 



Together with Eq. ©, they combine to yield the follow- 
ing result for the dynamical spin susceptibility 



X(r) 



2-e Tre 

tan — 

47T 4 



1 



1 



(Kog) 2 



(il) 



The exponent agrees with that of Ref s. lalalHJ, llll \12,\m 

If Jk is strictly 0, Gb{u) vanishes. For small but finite 
Jk, Eq. (0, which (more precisely, the term in between 
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FIG. 2: Static and dynamical local susceptibilities for e = 0.1. 
The red solid line corresponds to the quantum critical point, 
9 = 9o- 



the two equalities) comes with the assumption f3± < 1, 
cannot be satisfied. A solution does however exist for 
/?i > 1, in which case Gg(w) no longer diverges. It now 
follows from G~j^{uj) ~ const + (xA _1 that 

01 = 1 + e/2. (12) 

The non-divergence of G^ reflects the irrelevant nature 
of the Kondo interaction at the bosonic fixed point. 

Critical fixed point: The Kondo and bosonic terms are 
of the same order when f3\ = 1 — e + a,\. In this case, 
matching the dynamical parts leads to 



Oil 



e/2; /?! = !- e/2; 



A? = 



2-(l + «)e 7T6 



4tt(K <? c ) 2 



tan — ; B\ = — 

4 ' 47riVoAi 



(13) 
tan— .(14) 



The critical bosonic coupling g c needs to be numerically 
determined by matching the static part: 1/ Jk — ^b (0) = 
0. The resulting dynamical spin susceptibility is: 



X(r) 



(1 + K,)e 7T£ 

tan — 

47T 4 



1 



1 



(K g c y 



(15) 



We make two observations at this point. First, for 
infinitesimal e, our result agrees not only with the re- 
sult of the e-expansion for the large-N model [c.f. the 
subsection containing Eq. J3J|] but also that of the e- 
expansion for the the physically relevant SU(2) single- 
channel model|ll|,|l2|]. This implies that the large-N limit 
captures the quantum critical behavior of the physical 
systems, even though it yields a multichannel behavior 
on the Kondo side. This is reminiscent of the effects of 
spin anisotropies in the single-channel models: though 
the bosonic fixed point with Ising anisotropy (whose lo- 
cal susceptibility contains a finite Curie constant |n|) is 
very different from its counterpart with SU(2) symmetry, 
77 is the same for the QCPs of the two cases[ll|. Second, 
Eq. iTT5|) goes beyond the e-expansion result described 
earlier in the sense that it is valid beyond infinitesimal e. 

We have also determined the exponents of the sub- 
leading terms for Gf and Gb , finding a-i = e and fo = 1 . 

We next consider approaching the QCP from the 
Kondo side. Defining T* to be the crossover scale 
where the spectral function Af(uj) changes from Af(u>) ~ 



FIG. 3: Static and dynamical susceptibilities for e = 0.9. 



---^ for M « T* to A f (u) ~ TJI ^ for \w\ » T*, 
and inserting the ansatz into I/Jk = £.b(c<j = 0) = 
So dTg {T)G f (T), we find T* - (g c - g) 2 ^. Correspond- 
ingly, the susceptibility diverges as g approaches g c : 

X {T,u = 0, 5 < g c ) ~ {g c -g)-\ 7 = 2(1 - c)/c. (16) 

Our results so far represent an essentially complete an- 
alytical solution at zero temperature. The solution at fi- 
nite temperatures is considerably more difficult to obtain. 
In the following, we resort to numerical studies. 

Numerical results: For low-energy analysis, we find 
it important to solve the integral eqs. I|8I9|I in real fre- 
quencies and using a logarithmic discretization [supple- 
mented by linearly-spaced points; typically 250 points 
in total for the frequency range (—10,10)]. The con- 
vergence is determined in terms of j4b(w) and Af(uj): 
for each, the sum over the entire frequency range of 
the relative difference between two consecutive iterations 
is less than 10~ 5 . We choose k = 1/2 and No(u>) = 
(l/7r)exp(— ui 2 /ir). The nominal bare Kondo scale is 
T% = (l/JV )exp(-l/Ar J K ) w 0.06, for fixed J K = 0.8. 
The bosonic bath spectral function [cf. Eq. J2J] is cut off 
smoothly at A w 0.05. K% is fixed at T(2 - e). 

In Fig. |2Ji), we show the static local spin susceptibil- 
ity x(T) as a function of temperature for e = 0.1. For 
vanishing and relatively small bosonic coupling 3, x{T) 
saturates to a finite value in the zero-tempcrature limit 
- which characterizes the Kondo phase. At a critical 
coupling|22| of g/Tft « 56, x(T) diverges in a power-law 
fashion. Beyond this threshold coupling, x(T) continues 
to have a power-law divergence. Its amplitude, on the 
other hand, decreases with increasing g. These features 
are not restricted to small e; similar results are given in 
Fig. 3a), for e = 0.9. They are consistent with the analyt- 
ical results: Eq. (|16|) states that x(T = 0) should diverge 
as g approaches g c from below; for g > g c , on the other 
hand, Eq. l(TT|) implies that the susceptibility amplitude 
should indeed decrease as g increases. 

The quantum phase transition can also be seen in the 
dynamics. In Fig. 2b), we show the imaginary part of the 
dynamical local spin susceptibility, x"( w )j at the lowest 
studied temperature (T = 1Q- J T%) and for e = 0.1. On 
the Kondo side, it vanishes in the zero- frequency limit 
with an exponent close to 0.33, consistent with the ana- 
lytical result 2ai — 1 = 1/3. At g Cl it displays a power-law 
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FIG. 4: uj/T scaling of the dynamical spin susceptibility at 
g = g c , for e = 0.9 and k = 1/2. The scaling exponent is 0.1 



with Ising anisotropyll3|. Our results are also relevant 
in two other contexts. First, some generalized version of 
the Bose-Fermi Kondo model may capture the physics of 
impurities in high T c superconductors|2^ |25j. Results 
similar to what we report here should shed light on the 
dynamical properties of such systems, which have not yet 
been systematically addressed. Second, the competition 
between Kondo and paramagnon couplings also appears 
in the context of disorder in nearly magnetic metals |2fij. 
We acknowledge the support of NSF Grant DMR- 
0424125 and the Robert A. Welch Foundation (LZ, SK, 
and QS), DFG (SK), CNRS and Ecolc Polytechnique 
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well as the Lorcntz Center, Leiden and KITP, UCSB. 



divergence, with an exponent close to 0.9. On the bosonic 
side, the exponent remains the same while the amplitude 
decreases as g increases. Both exponents agree with the 
analytical result of 1 — 2a\ = 1 — e = 0.9. Similar results 
for e = 0.9 are shown in Fig. 3b). 

We now address the finite frequency and temperature 
behavior at the QCP in some detail. Fig. 21 plots the 
dynamical spin susceptibility for e = 0.9, clearly mani- 
festing an ui/T scaling. Remarkably, the scaling covers 
an overall 15 decades of lu/T (from IO -8 to about IO 7 )! 
For each temperature, the result falls on the scaling curve 
until an u of the order T^-. Similar behavior is observed 
for other values of e in the range < e < 1. The uj/T 
scali ng r eflects the interacting nature of the critical fixed 
point J23J: there is no energy scale other than T; the re- 
laxation rate, originating from some relevant coupling, 
has to be linear in T. 

The small but visible deviation from a complete col- 
lapse in the range lo,T < T^ reflects the influence of 
subleading contributions. For e = 0.9, a similar degree of 
uj/T scaling is also seen in the B- and /—spectral func- 
tions Ab{w, T) and Af(u>, T). The imaginary part of the 
/-self- energy, E'i(w,T), on the other hand, starts to de- 
viate from ui/T scaling at a smaller frequency, reflecting 
a stronger effect of the subleading terms. For e = 0.1, 
Ab{oj,T) and H"Auj,T) are less influenced by the sub- 
leading terms than Af{u),T) and x"(w,T). 

Our work serves as a basis for studies of the spin- 
isotropic Kondo lattice systems: by establishing the 
first non-perturbative approach to the Kondo-destroying 
QCP of the Bose-Fermi Kondo model, we are in the posi- 
tion to study the lattice problem in as systematic a way 
as the Quantum Monte-Carlo study of the Kondo lattice 
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